NASA Technical Memorandum 107538 



Scaling Methods for Simulating Aircraft 
In-Flight Icing Encounters 

David N. Anderson 

Lewis Research Center, Cleveland, Ohio 


and 

Gary A. Ruff 

Drexel University, Philadelphia, Pennsylvania 


Prepared for the 

Second International Symposium on Scale Modeling 
sponsored by the International Scale Modeling Committee 
Lexington, Kentucky, June 23-27, 1997 


National Aeronautics and 
Space Administration 


Lewis Research Center 


October 1997 


Available from 


NASA Center for Aerospace Information 
800 Elkridge Landing Road 
Lynthieum, MD 21090-2934 
Price Code: A03 


National Technical Information Service 
5287 Port Royal Road 
Springfield, VA 22100 
Price Code: A03 


SCALING METHODS FOR SIMULATING AIRCRAFT IN-FLIGHT ICING 
ENCOUNTERS 

David N. Anderson 

NASA Lewis Research Center 

Cleveland OH 44135 

and 

Gary A. Ruff* 

Mechanical Engineering and Mechanics Department 
Drexel University 
Philadelphia PA 19104 


ABSTRACT 

This paper discusses scaling methods which permit the use of subscale models in 
icing wind tunnels to simulate natural flight in icing. Natural icing conditions exist 
when air temperatures are below freezing but cloud water droplets are super-cooled 
liquid. Aircraft flying through such clouds are susceptible to the accretion of ice on 
the leading edges of unprotected components such as wings, tailplane and engine 
inlets. To establish the aerodynamic penalties of such ice accretion and to determine 
what parts need to be protected from ice accretion (by heating, for example), 
extensive flight and wind-tunnel testing is necessary for new aircraft and 
components. Testing in icing tunnels is less expensive than flight testing, is safer, 
and permits better control of the test conditions. However, because of limitations on 
both model size and operating conditions in wind tunnels, it is often necessary to 
perform tests with either size or test conditions scaled. This paper describes the 
theoretical background to the development of icing scaling methods, discusses four 
methods, and presents results of tests to validate them. 


NOMENCLATURE 

A c Accumulation parameter, dimensionless 

b Relative heat factor, dimensionless 

c Characteristic model length, cm 

c PtW Specific heat of water, cal/g K 

d Water droplet median volume diameter, pm 

Ko Modified droplet inertia parameter of Langmuir and Blodgett, dimensionless 
LWC Liquid-water content of cloud, g/m 3 
M Mach number, dimensionless 

n Freezing fraction of impinging water, dimensionless 

p Ambient static pressure, nt/m 2 


♦Summer Faculty Fellow at Lewis Research Center. 
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Re Reynolds number, dimensionless 

t Ambient static temperature, °C 

V Free-stream airspeed, m/s 

We Weber number, dimensionless 

<p Water-energy transfer term, °C 

Af Latent heat of freezing for water, cal/g 

6 Air-energy transfer term, °C 

Pi Density of ice, g/m 

Subscripts 

R Reference or full size 

S Scale 


INTRODUCTION 

Aircraft are susceptible to ice formations on engine inlets, tail planes and wings whenever 
flight through clouds at below-freezing temperatures occurs. Suspended water droplets in 
clouds are frequently super cooled; that is, the water exists as a liquid at a temperature 
below freezing. Supercooled water striking aircraft surfaces freezes, and the resulting ice 
accretions can have a significant and dangerous effect on the aerodynamic performance of 
an aircraft. In particular, ice formations decrease the lift and increase the drag. Large 
transport aircraft are protected against ice on critical components by directing hot air bled 
from the engine compressor to keep those surfaces warm enough to vaporize water. 
Some less critical surfaces may not be protected, and smaller aircraft often use 
intermittent impulse methods to remove small amounts of ice repeatedly. 

Aircraft and component manufacturers must thoroughly test new products to determine 
the effect of icing on their performance. This testing is performed both during the design 
process and for certification purposes. Flight testing is ultimately required but is 
expensive and can only be done when atmospheric icing conditions exist. Icing wind 
tunnels can simulate natural icing with water-spray and refrigeration systems and provide 
control of cloud conditions, temperature and airspeed to permit safe, convenient and 
relatively inexpensive testing. Some measurement of lift and drag changes can be made 
in the icing tunnel, and ice shapes are often recorded. More precise aerodynamic-penalty 
studies can be made in flight or in an aerodynamic tunnel by attaching a wood or 
styrofoam reproduction of the ice shape to the leading edge of the airfoil. 

Because of size limitations some components cannot be tested full size in an icing wind 
tunnel; furthermore, every tunnel has some bounds on the ranges of test conditions 
available for testing. For these reasons, it is desirable to establish reliable methods to 
scale either model size or test conditions. Efforts to establish scaling methods for icing 
tests began in the 1950’s and continue to the present. In general, to test scaling methods 
an ice shape is recorded for a reference, or full-size, condition, the scaling equations are 
applied to find the appropriate scaled condition, and the scaled ice shape is recorded. The 
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two ice shapes are then comparech If size has 

by multiplying the ctjonlmaKS of t^e scal^^P ^ ^ on subjective judgment. 

Whether the two shapes are in agreem q nerfect has always been difficult to 

and the quality of agreement when the match ' ' ' ^Ichve ice-shape comparisons, 
define, in this paper, in addition to the “"'^ “/“cf sL^ whtch pe Jts more 
we also present a new a^proac o qu y g lh<;oretical basis for traditional 

^mXT^ss bTefly a Buckingham-rt analysis, descnbe four scaitng methods 
S present results from a series of tests of these scaitng methods. 


DERIVATION OF SCALING EQUATIONS 


The traditional approach to 

similitude in the geometry, the flowfi , J I „ have been denV ed 

catch, and the energy balance at the surfi^V ^ sc^ g ^ ^ ^ 

,1.23,4,51 which provide s—^-me^or water feezes on 

accretion. When ice accret , surface energy balance have been shown 

impact, simple scaling metho s w ic lg ^ accredon f or which water does not 

to work successfully [6], On the other , \ c ^ al&d by methods which include the energy 
freeze immediately on impact, can y ohvsics of ice accretion have 

balance. Recently, attempts to understand *£££, effect 0 „ the final 

shown that for glaze ice, surface p ch f the similitude requirements which 

ice shape [7,8]. In this section, we will address eacn oi 

have been satisfied by traditional methods. 

Geometry The scaled and reference bTI 

entire surfaces. It is assumed that as 1C ® ’ An alternat e approach is currently being 

scaled representation of the reference ices P _ geometry is full-size while the 

studied for airfoils in which t e ea lI \ 8 Th design of the truncated airfoil is 

remainder of the seale airfoil ,s truncated 9 'Ol jbe destgn me(hod wil , 

such that the flowfield around it matches that of the reference a, no 

not be discussed in this paper. 

Flowfield The flowfield over the 

OTer ‘ kas ‘ u,e lead,ng - ed8e region of 

the scale model match those for the reference. 

Droplet Trajectory To insure that the massif both 

^r.^1 water catch have to he similar. The water catch will 
be discussed in the next section. 
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[1 1J. PUbHshed in Ref ‘ W ™ d 

reference size, d R , according to Eq. (1): r ° P Ct S1ZC ’ ds ’ muSt re,ate to the 
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^ o“ t z;ri: ( z 7oir,T ly rr for the range ° f *«*« 

effect on the drag of a moving soherLl i ■ T t""* 1 by a W> ro ™i"S the Re 
somewhat more accurate approach was used by Rufn’jP h d*, “ ex P ression A 

size which satisfied equality of the modified in „ 5 ’ he de,ermm ed the scale droplet 

add Blodgett |12] defined the Id^d fnLt par ™ eKrs ' «u = Langmuir 
inertia of a spherical droplet Either En rn p “' ra ''' cter 10 lnc l“de effects °f both Re and 
for compu ^calculation “ *» -»y ptogrammed 

fecz^T eventually C “^“on" “» «*. surface is assumed ,0 

dimensional accuAlat.onT™, Z^WCV rT “"J*. d ™ * 

Will accrete the same relative quantity of ice the - 1 ° P " T°.' nSUre that the sca,e test 

to the reference. Thus, ’ SCa 6 accumu| ation parameter is matched 


— = f— y 

t r l C R A 


V 


LWC ’ 


LWC s j{V s 


( 2 ) 


f r s T ediateiy ° n impact ’ * is ° ni ^ ——y 

sin.ili.ude in energy ^ 

and freezing was performed by^Snger H 3^a^d tl" b^^lT’h impingement 

methods since that time. Messineer 1,16 basis of most scaling 

surface due to convection, ice sublimation wate 3 " 06 mC uded the ,oss of heat from the 
required to warm impinging water to the VreezTnrrem^rato^The a " d hea ‘ 

surface is due to release of latent heat on feezing aTS, 7 A 8 ° f heat aI the 
water droplets. Ruff [5] added terms for th» a 8 d the kinetic energy of incoming 
and for hea, earned ° f ^ thro “ 8h “« -*.* 

balanced the^fr^zing' fraction, “» " energy 

impingement region; the relative heat factor b^isth ° f Wa j. er wh,ch freezes within an 

impinging water to the ability of the airflow tn * u ° f ^ t0tal heat ca P acit y of the 
parameters used are 0 wWch is a ^ ^ 1,16 SUrface ‘ Tw « other 

transfer, and ft which groups the air-energy Lsfer^TIh’'" 1 ' dr0P ‘ e ‘ 
written using these parameters, it becomes WhCn he Cnergy balance is 
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p,w 
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Various scaling methods have selected one or more of the parameters n, <p, 6 , and b to 
match between scale and reference values. 


BUCKINGHAM-n ANALYSIS 

The scaling parameters discussed in the preceding section result from a 
phenomenological approach in that they are derived from a set of equations that describe 
our current understanding of the ice-accretion process. Classical scaling by dimensional 
analysis and application of the Buckingham-7t methodology has also been applied to icing 
scaling. The premise behind this approach is that if the proper dimensionless parameters, 
or 7t terms, can be identified, any one of the 7t terms can be written as a function of the 
remaining 7t terms. Only the parameters relevant to physical phenomena and their 
dimensions must be specified. This approach is very attractive because it requires 
minimal knowledge about the physics of icing, and we presently lack a full understanding 
of all the physical processes involved in ice accretion. 

This methodology was applied to icing scaling by Bilanin [7] who identified 23 variables 
that play a role in the ice-accretion process. With 23 variables and 4 dimensions (length, 
mass, time and temperature), there are 19 possible dimensionless parameters. Three of 
the parameters are easily recognizable as the Mach, Reynolds and Weber numbers. 
Others are related to the drop size, inter-droplet spacing and ffee-stream temperature. 
Still others are ratios of the physical properties of water and ice which will, of course, be 
matched automatically between scale and reference situations. It’s possible to show that 
all of the traditional icing similitude parameters discussed previously can be expressed as 
functions of some of the 7t parameters. If all 19 of the re parameters for the scale test were 
simultaneously matched to their respective reference, or full-size, values the ice accretion 
should be rigorously scaled. 

Unfortunately, not all n parameters can be simultaneously matched from scale to 
reference conditions. For example, a constant M requires approximately that V s = V R , a 
constant Re requires that c s /c R = V R /V S and a constant We requires that c s /c R = (V R /V S ) . 
Except for the special case when c s = c R , these are inconsistent restrictions. Fortunately, 
it may not be necessary to match all 7i parameters to scale successfully; the methods 
described in Refs. [1] - [6] were based on the traditional icing similitude parameters 
which ignored many of the 7t parameters, yet scaled the geometry of ice accretions fairly 
accurately for some conditions. Although Bilanin’ s Buckingham-Jt analysis has not 
proved to be practical as a scaling method by itself, it has served to identify parameters 
which may have been overlooked by the traditional methods. One such parameter is the 
We which has now been incorporated into a new scaling method discussed below. 
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SCALING METHODS 


Icing scaling methods have matched the scale and reference values of a variety of the 
terms discussed above to find the scale test conditions. The scaling methods will be 
described next in terms of the parameters matched. A summary of the 4 scaling methods 
and the similitude parameters which they satisfy is given in Table 1. This list of scaling 


Table 1. Similitude Terms Satisfied by Four Scaling Methods. 


Method 

Drop 

Traj 

Drop 

Catch 

Rel. 

Heat 

Fact 

Freez. 

Fract. 

Drop 

Engy 

Trans 

Air 

Engy 

Trans 

M 

Re 

We 

LWC x time 

X 

X 



X 

X 

X 

X 

X 

Olsen 

X 

X 


X 



X 

X 

X 

Ruff (mod) 

X 

X 


X 

X 





Const- We 

X 

X 


X 

X 




X 


methods is not complete; in particular, the French method [3] will not be discussed here. 
It was evaluated in Ref. [6] and has been used extensively in the ONERA icing tunnel at 
Modane. 

“ LWC x time = constant” This is the simplest and oldest scaling “law.” It applies 
when the model size is not scaled and all parameters except the liquid-water content can 
be matched to the desired test conditions. In this situation, the law states that the amount 
of accreted ice for the scaled LWC will be the same as for the desired LWC if the product 
of accretion time and LWC for the scaled test equals that for the desired, or reference, 
encounter being simulated. This expression is derived from Eq. (2) with c s = c R and V s = 
Vr. In addition, with ps = p R , Eq. (1) gives ds = d R . This method requires that the static 
temperatures also match; i.e., = t R . As a result, and because of the other matched 

parameters, tps = <f> R , 6 S = d R , Wes = We R , Re s - Re R and Ms = M R as well. 


Olsen A variation on “ LWC x time = constant” is the Olsen method [6, 14] in which, 
again, only the LWC is to be scaled. As in the case of the “LWC x time = constant” 
approach, matching the model size and test conditions other than LWC results in a match 
of the We, Re and M as well. In this method, however, the scale static temperature does 
not equal the reference temperature, but is found by matching the scale and reference 
freezing fraction, n. 
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Ruff The last two methods to be discussed are used when the scale model size is 
different from the reference. The Ruff method was developed at the AEDC by Ruff 
[5] and is sometimes known as the AEDC method. It was intended originally for use 
in wind tunnels with altitude capability. The user first selects a scale model size and 
test velocity. The static temperature for the scale test can then be found by matching 
the water-energy transport term, (f> . The scale droplet size is determined by 
matching the modified inertia parameter, K 0 , which insures that the scale droplet 
trajectories will be the same as the reference case. The scale static pressure is found 
by matching the air-energy transport terms, 0. For tunnels which cannot control the 
static pressure independently of velocity, a modified Ruff method has been used in 
which the 6 terms are not matched. Next, the freezing fraction is matched to 
establish the scale liquid-water content, and, finally, the icing encounter time is 
found from matching the scale and reference accumulation parameters. 

Constant- We The Buckingham-7t analysis of Bilanin [7] showed that the Weber 
number may be a significant icing scaling parameter. This method assumes that We 
is of greater importance than either Re or M, as no attempt is made to match these 
two parameters to the reference values. The user chooses the model size; all other 
scale parameters are determined by similitude requirements. The 3 equations formed 
by matching the Weber number, the modified inertia parameter, and the water- 
energy transfer term are solved simultaneously to give the scale velocity, droplet size 
and static temperature. The scale liquid- water content and icing time are then found 
from matching the freezing fraction and accumulation parameter, respectively. 


EXPERIMENTAL METHODS 


Test Facility and Models Tests to validate scaling methods were performed in 
the NASA Lewis Icing Research Tunnel (IRT) shown in Fig. 1. It is a closed-loop 
tunnel with a test section 1 .8 m high by 2.7 m wide. Temperature can be controlled 
from -30°Cto 4°C. The water spray system gives a range of liquid-water content 

and water droplet 
size which covers a 
significant portion 
of the FAA Part 25 
Appendix C icing 
envelope. Veloci- 
ties of up to 160 
m/s are possible. 



Fig. 2 is a photo of 
a 53.3-cm-chord, 
1.8-m-span NACA 
0012 airfoil 

mounted vertically 
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Figure 2. 53.3-cm-Chord NACA 0012 Airfoil 
installed in IRT. 


in the IRT test section for tests to 
be reported here. This model had a 
uniform chord over the full span 
and was unswept. Test conditions 
were selected to represent 
reference cases, and the various 
scaling methods applied to 
determine the corresponding scale 
test conditions. When the scale 
tests involved a size change, an 
NACA 0012 model with 26.7-cm 
chord was used. Tests were run at 
both reference and scale 
conditions, two-dimensional cuts 
through the resulting ice accretions 
were made at the center of the 
tunnel test section, and ice shapes 
were recorded by tracing the ice 
outline onto a cardboard template. 
These tracings were then digitized 
for computer storage. From these 
computer files, the shapes were 
analyzed and compared. 


Quantitative Analysis of Ice-Shape Features In general, evaluations of icing 
scaling methods have relied on qualitative comparisons of the scale and reference ice 
shapes. With this approach, the experience, judgement and objectivity of the 
researcher determine to some extent whether a scaling method is acceptable. In 
practice, what constitutes an acceptably scaled ice accretion actually depends on the 
purpose of the test. In aerodynamic tests, it is most critical that the sub-scale 
accretion have the same lift and drag coefficients as the full-scale. In other 
applications, geometric parameters such as the width or mass of the accretion is most 
important. In still others, the critical parameter may be the mass of ice shed. For 
general evaluation of scaling methods, however, a comparison of ice shapes is most 

appropriate. In this study, 
both qualitative and 
quantitative comparisons of 
ice shapes will be made. 

Several characteristic 

dimensions representative of 
the overall shape of a typical 
glaze ice accretion were 
identified. As shown in Fig. 
3, these dimensions included 
the thickness of the ice ac- 
cretion at the stagnation 



Figure 3. Ice-Shape Characteristic Dimensions. 
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point the maximum thickness, the maximum width of the ice accretion the 
impingement width, horn length and horn angle. These characteristics wen 
measured on the main ice shape. Downstream of the primary glaze ice shape 
region in which rime feathers form. The features of these feathers vary considerably 
from test to test and from one location to another on the model. Because of is 
variability, the feather region is ignored in comparing ice shapes. Measurements o 
the characteristic dimensions were made by hand for this study. 

Ice-Shape Repeatability To be judged acceptable, scaling methods must 
produce ice shapes that are similar to the reference shape within the typical shape 
variability from run to run. To establish this variability, several test conditions w 
repeated and the shapes compared. In Fig. 4 (a) are shown results for tests made m 
October 1995, December, 1995 and June, 1996 at the same tunnel conditions. In 




(a) f, -12°C; V, 67 m/s; d, 30 pm; 
LWC, 1 g/m 3 ; T, 7.3 min. 


(b) t , -9°C; V, 89 m/s; d, 40 pm; 
LWC , .55 g/m 3 ; T, 10.0 min. 


Figure 4 . Repeatability of Ice Shapes. 53.3-cm-Chord NACA 0012 Airfoil. 

Fig 4 (b) are results from tests in October, 1995 and December, 1995 at another set 
of conditions The reference shape (October, 1995 test) is shown with a solid line in 
both parts of the figure. Although small differences in ice shape are apparent from 
this quahtative comparison, the 1RT generally gives fairly repeatable ice shapes. 

In addition to this subjective evaluation, a quantitative assessment of ice-shape 
in aaaiuun iu i repeatability was also made. 

. . The six characteristic 

Table 2. Variability of Six Characteristic dimensio ns of the ice-shapes in 

Ice-Shape Dimensions. pig 4 were measured and 

averages for C3.ch dimension 
obtained separately for Fig. 4 
(a) and Fig. 4 (b). The 
difference between each 
dimension and the average 
value of that dimension was 
then obtained. Finally, the 
absolute values of these 
differences were averaged and 


Ice Feature 

Average Percent Difference 
from Mean Dimension 
Fie. 4 (al Fig. 4 (b) 

Stag. Thickness 

8.7 

11.4 

Max. Thickness 

2.8 

5.2 

Max. Width 

10.4 

6.6 

Horn Length 

0.4 

4.2 

Horn Angle 

7.7 

23.1 

Impinge. Width 

9.4 

13.7 
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reported as a percent of the mean dimension in Table 2. For these two icing 
conditions, the maximum thickness and horn length had the best repeatability while 
the stagnation-zone thickness, maximum ice-shape width, horn angle and 
impingement width were less repeatable. Unfortunately, some judgement is still 
required to determine some of these measurements. The impingement width, which 
s dependent on the droplet trajectory, is often particularly difficult to define and 
this uncertainty was reflected m the relatively poor repeatability of this dimension. 

both 1113 6 ^fi ° f d ™ plet trajectories and impingement limits has been verified 
computationally and experimentally using temperatures above freezing [1 1,17], 

' , '! in r “ S ° n “ ble t0 “P" th « the characteristic dimensions of an ice shape can be 
defined better as more repeat data are available. In Table 2, for most dimensions, the 
average deviation from the mean dimension was less when based on three ice shapes 

Fk 4 fhvl S Th7 9 ! § ' 4 ? )} than WhCn ° nly tW ° ice Shapes were used (results for 
g t ki 3b e 2 3 S ° indlcates that most ice-shape characteristic dimensions are 
repeatable to about ±10%, with horn angle somewhat more difficult to repeat 
Ability to reconstruct most ice dimensions to within ±10% is thus a reasonable goal 
for scaling methods with somewhat more relaxed expectations to reproduce horn 


EVALUATION OF THE SCALING METHODS 


Fig al 5 for and SCale ice shapes « —pared in 

ig. 5 for the LWC x time = constant rule (Fig. 5 (a)) and for the Olsen scaling 



(a) Scaling Using "LWC x time = Constant". (b) Scaling Using Olsen Method. 


Reference 

Scaled to .8 Ref. LWC 


LWCxtime 

Olsen 

Parameters for 

Both Methods 

t 

t 

V 

d 

LWC 

T 

°C 

°C 

m/s 

pm 

g/m 3 

min 

-12 

-12 

67 

30 

1.00 

7.3 

-12 

-10 

67 

30 

.80 

9.1 

: -12 

-15 

67 

30 

1.40 

5.2 


Figure 5. Scaling Liquid-Water Content. 53.3-cm-Chord NACA 0012 Airfoil. 

method (Fig. 5 (b)). In each case, the reference test (solid line) was performed with 
iqui ~ water content of 1 g/m , and scale tests were made with LWC's of .8 and 1 .4 
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Table 3. Quantitative Evaluation of “ LWC x time = Con- 
stant” and Olsen Scaling Methods. Reference 
LWC, 1.0 g/m 3 ; Scale LWC, .8 and 1.4 g/m 3 . 


Ice Feature 

Percent Difference from Reference 


LWC x Time 

Olsen 


.8 g/m 3 

1 .4 g/m 3 

.8 g/m 3 

1 .4 g/m 3 

Stag. Thickness 

10.0 

-20.0 

0.0 

5.0 

Max. Thickness 

5.7 

-17.1 

0.0 

-4.3 

Max. Width 

-17.0 

30.1 

11.3 

-1.9 

Horn Length 

0.0 

4.6 

-2.8 

-0.9 

Horn Angle 

-43.5 

39.1 

-4.4 

0.0 


g/m 3 . In this test series, 
both methods of scaling 
produced an ice shape 
which simulated the 
reference shape fairly 
well when the liquid- 
water content was 
reduced to .8 g/m 3 , but 
the Olsen method gave 
a superior match when 
the scale LWC was 
increased to 1 .4 g/m . 


The ice shapes shown in Fig. 5 were analyzed quantitatively, and the results are 
reported in Table 3. Because of the difficulty in defining the impingement width, 
this dimension will not be included. With the exception of the horn angle, the scaled 
dimensions resulting from using the “LWC x time = constant” method were close to 
being within the acceptable limit of ±10% of the reference dimensions when the 
liquid-water content was scaled from 1 to .8 g/m 3 . However, for the 1.4-g/m 
scaling case, this method produced an ice shape with dimensions significantly 
different from the reference. 


The Olsen method provided scaled shapes whose dimensions closely matched the 
reference for both scaling cases. The match of the horn angle produced by the Olsen 
method is particularly notable. The formation of horns depends on the dynamics of 
liquid water on the surface of the ice accretion. The success of the Olsen method 
over the “LWC x time = constant” method indicates that the freezing fraction, 
matched to the reference value by the Olsen method, is of greater importance in 
determining final ice shape than the air-energy or water-energy transport terms, 
which are matched by the “LWC x time = constant” method. 

Scaling Size Results of tests with size scaled to Vz the reference value are shown in 
Fig. 6. The solid line in each part of the figure represents the reference ice shape 
which was the same for each scaled test. The dotted-line shape is the scaled result. 
Both the Ruff and the constant- We methods gave liquid- water contents for the Vz- 
scale conditions which were outside the range of the operating map for the IRT; 
consequently, the LWC for the scale tests was selected to be 0.8 g/m 3 , and Olsen 
scaling was applied to maintain the same freezing fraction as the reference test. 
Thus, both size scaling and LWC scaling were combined for this evaluation. 

Because the IRT does not provide control over the test-section pressure, the 
modified form of the Ruff method was used. The results in Fig. 6 (a) show that the 
scale ice shape was of similar size to the reference although, qualitatively, the 
stagnation thickness and horn location were slightly different. The ice tracings 
shown in Fig. 6 (b) for the constant- We scaling conditions appear to match the 
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(a) Scaling Using Ruff (Mod) Method. 

c 



cm 

- Reference 

53.3 

• Ruff Method 

26.7 

• Constant- We Method 

26.7 


(b) Scaling Using Constant- We Method. 


t 

V 

d 

LWC 

r 

°C 

m/s 

pm 

g/m 3 

min 

-12 

67 

30 

1.00 

7.3 

-8 

67 

20 

.80 

4.6 

-10 

88 

18 

.80 

3.5 


Figure 6. Model Size Scaled to Vi Original Chord. NACA 0012 Airfoils. 


Table 4. Quantitative Evaluation of Modified 
Ruff and Constant- We Methods for 
1/2-Size Scaling. Test Conditions 
Given in Fig. 6. 


reference shape slightly better. 
The quantitative comparison of 
the dimensions shown in Table 4 
show that the constant- We method 
provided agreement with the 
reference shape which was only 
slightly better than that for the 
modified Ruff method for this 
case. Because the constant- We 
method places more restrictions 
on the scaling conditions, it’s not 
surprising that it was fairly 
successful; however, the fact that 
the Ruff method, which ignores 
We, was relatively successful 
suggests that We may not be of 
major importance to icing physics. Clearly, there is opportunity for additional 
improvement in icing scaling by better understanding the phenomena involved in the 
ice-accretion process. 


Ice Feature 

Percent Difference from 
Reference 

Ruff (Mod) Const- We 

Stag. 

Thickness 

25.0 

25.0 

Max. 

Thickness 

-17.1 

-5.7 

Max. Width 

-22.6 

-18.9 

Horn Length 

-14.8 

-13.0 

Horn Angle 

52.2 

39.1 


CONCLUDING REMARKS 

This paper described the theoretical background leading to the development of four 
icing scaling methods. The phenomenological basis for current scaling technology 
was presented and compared to the classical approach using the Buckingham-^ 
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methodology. Tests to evaluate these scaling methods were conducted in the Icing 
Research Tunnel at NASA Lewis Research Center, and results were presented. The 
“LWC x time = constant” and Olsen methods can be used to scale test conditions 
while the constant- We and modified Ruff methods scale test article size. A method 
to quantify the goodness of the scale ice shapes based on measuring six 
characteristic dimensions was proposed. Results from this quantitative approach 
supplemented visual comparisons performed by overlaying two-dimensional tracings 
of the ice accretion. The conclusions from this study are summarized below: 

1. Quantitative verification of icing scaling methods is helpful in defining their 
accuracy and increasing confidence in their use. 

2. When comparing the characteristic ice-shape dimensions from repeat conditions, 
differences on the order of ±10% can result if only two tests are compared. 
Characteristics such as horn angle and impingement width are sometimes 
difficult to define, resulting in large run-to-run deviations in these dimensions. 
Improved definition of dimensions resulted from a greater number of samples. 

3. When liquid-water content was scaled using a full-size model, the Olsen method 
produced better results than the “LWC x time = constant” method. The former 
maintains the freezing fraction constant between reference and scaled conditions, 
while the latter matches the water-energy-transfer and the air-energy-transfer 
terms. This result suggests that the freezing fraction has a greater effect on ice 
shape than the water- or air-energy-transfer terms. 

4. For scaling to x /i size at the conditions tested, the constant- We method produced a 
slightly better match of ice shape than the modified Ruff method. Either of these 
methods appeared to provide at least approximate scaling for ‘/ 2 -size models. 
However, the relative success of the less-restrictive Ruff method suggests that 
We may not be as important to ice-accretion physics as once thought. Additional 
study of the parameters of most importance to the development of ice shapes is 
needed to improve scaling methods further. 

Several directions for future work are evident from these results. The quantitative 
evaluation of ice accretion shapes appears to be a promising tool, but the observed 
variation between any two repeat conditions showed that a large number of repeated 
tests providing better statistics are required to perform more detailed evaluations of 
icing scaling methods. Characteristics of the ice shape other than those discussed 
here may also need to be considered for future quantitative evaluations. To improve 
scaling methods, other parameters, such as Re, need to be examined to determine 
their importance relative to We. Finally, parameters relevant to the dynamics of 
liquid water on the surface of an ice accretion during the freezing process should be 
identified and included in icing scaling methods. 
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